Spin-flux phase in the Kondo lattice model with classical localized spins 
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We provide numerical evidence that a spin-flux phase exists as a ground state of Kondo lattice 
model with classical local spins on a square lattice. This state manifests itself as a double-Q magnetic 
order in the classical spins with spin density at both (0, 7r) and (n, 0) and further exhibits fermionic 
spin currents around an elementary plaquette of the square lattice. We examine the spin-wave 
spectrum of this phase. We further discus an extension to a face centered cubic (FCC) lattice where 
a spin-flux phase may also exist. On the FCC lattice the spin-flux phase manifests itself as a triple-Q 
magnetically ordered state and may exist in 7-Mn alloys. 



The Kondo lattice model with classical local spins has emerged as one of the simplest models that can account 
for some of the physics of the manganites |jj and the cuprates For the manganites the ferromagnetic Kondo 
lattice model gives rise to the double exchange model which has been argued to be the relevant model to explain the 
physics in these materials [pi. For the cuprates the assumption of classical local spins is clearly unrealistic, however 
the antifcrromagnetic Kondo lattice model gives rise to many insights into high T c materials. For example, it has 
been used to understand the appearance of d-wave superconductivity |^|-|(| and it also gives rise to incommensurate 
magnetic and stripe Q structures that have been experimentally observed ^,|) . One aspect of this model that has 
recently gained interest is the appearance of a Berry phase in the fcrmion wave function that arises when the fermion 
spin is strongly pinned to the local classical spin orientation []1 Of - This Berry phase has been argued to give rise to a 
flux phase ground state on a square lattice in the manganites to an anomalous hall-effect in ferromagnets [fL2[ , 
and to a quantized hall conductance in Kagome lattices E3). In this paper we examine the conditions under which 
this Berry phase gives rise to novel ground state structures. In particular, we give numerical evidence that a spin-flux 
phase appears as a ground state structure of this model on a square lattice. This phase is analogous to but quite 



different from the flux phases that are usually discussed in the context of the cuprates 14 15[| . The difference arises 
because the latter flux phases exhibit a finite current around each elementary plaquette of the square lattice, but 
in our case spin currents exist (for which the up and down electrons have opposite currents around an elementary 
plaquette). On the square lattice the phase discussed here has a spin-flux of it through each elementary plaquette. In 
this regard, it is of interest to note that a n spin-flux phase has been argued to be central in explaining the normal 
state properties of the cuprates by John and co-workers 

In this paper we will first demonstrate numerically that a double-Q magnetic structure exists as a ground state of 
the Kondo lattice model. We then demonstrate that such a state is a spin-flux state with circulating spin currents 
and estimate the stability region of this phase. We also determine the spin-wave spectrum of the double-Q magnetic 
phase and demonstrate that a spin flux phase may exist for this model on a face centered cubic (FCC) lattice. The 
spin- flux phase on the FCC lattice manifests itself as a triple-Q magnetically ordered state and has a spin flux of w/2 
through each elementary triangular plaquette that lies in the planes having Miller indices (1,1,1) (and equivalent 
symmetry planes). 



The model we study here is 



H = -* E ( C L C J« + -J^Si-Si + J'^SrSj, (1) 

where Cj a creates an electron at site i = (i x ,iy) with spin projection a, Si=^ Q/3 C; Q ,cr Q ,^Ci / 3 is the spin of the mobile 
electron, the Pauli matrices are denoted by <x, Si is the localized spin at site i, (ij) denotes nearest-neighbor (NN) 
lattice sites, t is the NN-hopping amplitude for the electrons, J is a coupling between the spins of the mobile and 
localized degrees of freedom, and J' > is a direct AF coupling between the localized classical spins. Throughout 
this article the unit of energy will correspond to t = 1. For the numerical studies a Monte Carlo technique was used. 
This involves no "sign problems" so that by this procedure temperatures as low as T=0.005 at any density can be 
reached. The present study has been performed mostly on 6x6 lattices with periodic boundary conditions (PBC), 
but occasional runs were made also using open and antiperiodic BC as well as different lattice sizes (up to 12x12 
lattices). The specific numerical technique used here involves a standard Metropolis algorithm for the classical spins 
and an exact diagonalization for the itinerant electrons. The details of the method are described in Ref. fll7||. 



1 



The spin-flux phase was identified numerically by studying the classical spin structure factor which is the Fourier 
transform of the static spin-spin correlation function S(q) = J2 n m e iq ^ n ~ m ^(S n • S m ). In particular, it was found 
that for various JS and J'S 2 in the vicinity of electron density (n) = 0.5 the structure factor was peaked at Q = (0, 7r) 
(Q y ) and Q = (71", 0) (Q x ) (see Fig. 0) To understand the possible ground states for a spin density with this wave 
vector it is useful to look at a Ginzburg Landau free energy. The order parameter is determined by the two vectors 
Mo iTr and M^o- The free energy can be simply constructed by noting that the relevant space group representation 
transforms as a vector under spin rotations and as a scalar (that is as an A\ g representation) under the little co-group 
D2H of the wavevector Q = (0,7r). The most general dimensionless Ginzburg Landau free energy is 

F = -(Ml„ + Ml fi ) + (M^ + M 2 ^) 2 + /?M^M 2 !0 + ft(M„ ; „ ■ (2) 

The minimization of this energy leads to three possible ground states: (a) (M 0iTr , M^.q) = (M, 0), (b) (Mq.^M^.q) = 
(Mi, M 2 ) with Mi -M 2 = 0, (c) (M 0)7r , M^q) = (M, M). The double-Q phase (b) we will argue below is the spin-flux 
phase which in fact corresponds to a particular representation of flux phase proposed by Yamanaka et al. JTT|] (note 
that this phase does not lead to a peak in S(q) at (w/2,tv/2) as suggested in Ref. 0). The double-Q phase (c) 
corresponds to ordering only one half of the local moments and is therefore not a likely ground state for this model 
(note however that there exists numerical evidence for this phase in a periodic Anderson model on a square lattice 
jl9| ) . To distinguish numerically which of these three phases corresponds to the phase found here the spin correlations 
were examined by evaluating Si-Sj = S 2 cos (the spin dot product of NN spins) for each pair of NN spins and 
plotting the value of cos 9^ in a histogram. The results are shown in Fig. |l| for JS = 2 and J'S 2 — 0. From this 
figure it is clear that NN spins are orthogonal which implies the double-Q order of phase (b) above. 
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FIG. 1. (a) Structure factor S(q), (b) histogram of snapshots for the nearest neighbor spin dot product, (c) and resulting 
classical spin structure for JS = 2, J' S 2 = 0, T = 1/200, and (n) = 0.5 on a 6 by 6 lattice with periodic boundary conditions. 

To understand the electronic properties of this double-Q magnetic phase we fix the classical spins and find the 
fermion eigenstates (for the double-Q state this is reasonable because the spin structure factor is very strongly peaked 
at (0, 7r) and (ir, 0) with little weight at other q values as can be seen in Fig |l|). The classical spin orientation is given 
by Si = (5/2)[(— l) 4x + (— l) ly , (— l) 4x — (— l)* a ,0]. Solving for the eigenstates of the resulting electronic Hamiltonian 
results in four bands with dispersions 



e k = ±\ (JS) 2 + 4(cos 2 k x + cos 2 k y ) ± 2J2(JS) 2 (cos 2 k x + cos 2 k y ) + 16 cos 2 k x cos 2 k y 



(3) 



where (k x ,k y ) £ { | k x + k y \ < tt} n { | k x — k y | < tt} are restricted to one half the original Brillouin zone. The density of 
states (DOS) is linear in |fe| for (n) = 0.5 which is characteristic of the Dirac spectrum that appears for 7r-flux phases 
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Jl4] , fil| . Also note that the dispersion relation is independent of the sign of J; consequently if this flux phase is the 
ground state for positive J then it must also be the ground state for negative J. In the limit J = oo the dispersion 
reduces to that found in Ref. O], To identify this phase as a spin- flux state the spin current from site i to j was 
determined 



it 



a,/3 



(4) 



It was found that only [jij] z is non-zero and it is non-zero only for NN sites. The resulting spin currents circulate 
neighboring plaquettes in opposite directions. The charge current was found to be zero. This spin current pattern 
implies a spin-flux of ir exists in each elementary plaquette. Note that in the limit J = oo this result is intuitively 
clear; in this limit the spin of the electron is tied to the local moment so that when the fermion travels around a 
plaquette the spin changes by 2tt which implies that the wavefunction changes sign. 

The phase diagram for (n) = 0.5 as a function of l/(JS) and J'S 2 is shown in Fig. |[ The solid phase boundaries 
were found by comparing the energy of the flux phase to that of the canted magnetic and spin density wave (SDW) 
phases (the energies of the helical SDW phases agree with those found in Ref. (20|). At larger J'S 2 the flux phase is 
found to be unstable to a SDW phase characterized by Si = Sv2(— l) 1 " cos(i x 7r/2 — 7r/4) (note this is not a helical 
SDW state). This agrees with the structure found numerically. 




FIG. 2. Phase diagram for (n) = 0.5. The filled squares represent the numerical results on a 6 x 6 lattice with T = 1/200 
and with periodic boundary conditions. 



It is of interest to determine the spin-wave spectrum arising from the spin-flux phase. This can be done by using 
the spin-wave approximation that was introduced by Kubo and Ohata [ pl| and later used by Furukawa for the double- 
exchange model |22| . The local spins are described by a local co-ordinate system in which each classical spin is aligned 
along the z direction and the spin- wave operators Sf ~ V2Sai, S^ ~ V2Sa\, and Sf = S — a\a,i are introduced. 
The spin- wave spectrum is found by keeping all 1/S corrections to the magnon self-energy. Here we consider the limit 
J = oo. This results in the following effective boson Hamiltonian for the spin waves 



^~][II(k)a£a k + A(k)a k a_ k + h. 



(•5) 



with k summed over the whole Brillouin zone of the square lattice, 



n(k) 



2SN ^ 
q 



cos( 



9q — #k+q)-Ek+q ~ 



E„ + i?kH 



[1 - cos(5#k+ q - q )] > + J'S(cosk x + cosky), (6) 



A(k) 



2SN 



-Efc+q-Eq 
Eq + -Ek+c 



[C0s(2$q + 2#k+q) — COs(# q — #k+q)] + J'S(cOSk x + COSfcy), 



(7) 



cos 9\, 



[cosk x + cosky]/ Ek, sin 9^ = [cosk x — cosky]/ E^, and 



2 J cos 2 k x + cos 2 k y . The spin 



wave dispersion is given by cj£ = -\/n(k) 2 — |A(k)| 2 . The eigenstate of this mode is given by 5S q (r) = 



e *Q- r (n(q)+A(q))^ 



e lC V r (n(q) + A(q))M (recall the ordered moment is S (r) = e^ r ^- 
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The resulting spin- wave spectrum (see Fig. 3) agrees with the general form required by phenomenological 



arguments (found by using the method of Zhu and Walker E3]). 
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FIG. 3. Magnon dispersion relation for the double-Q magnetic phase (J' S 2 = 0.1). 

Given that the spin-flux phase was found to be stable on a square lattice it is natural to ask whether such 
states can be realized on other lattice structures. We argue that a spin-flux phase is likely to be stable on 
a FCC lattice. For an FCC lattice two degenerate ground states of the NN classical Hciscnberg model arc 
S 3 q(R) = [(-l) h+h ,(-l) l2+ls ,(-l) h+l3 ]/V^ and S 1Q (R) = [(-l) /l+ ' 2 ,0,0] where the FCC lattice is spanned by 
R = [l%a(x + y)/2, l^a{x + z)/2, ha(y + z)/2] (note that there exists a continuous degeneracy in the ground state, but 
in the presence of the Kondo coupling only two states are relevant). In the limit J = oo the structure S3Q gives rise 
to a spin-flux phase with the spectrum 



4 / 2 k * 

£k = ± VlV cos T cos 2 



2 k v 



k h 
• 2 Z 9 V 

sm — cos^ — 
2 2 



• 2 n 

sm — cos^ — 



(8) 



where the momenta are restricted to the region of the Brillouin zone where — 2ir < k z < 0. Note that for 

fc. 



= 0.5 the 



DOS is again linear in energy. For the structure Siq the dispersion is = — 4£ cos cos Assuming that one of 
these ground states is stable (that is taking J' to be sufficiently large) then it is found that at (n) = 1 the Siq state 
is stable while at (n) = 0.5 the S3Q state is stable. There is a transition between these two states at (n) = 0.7. As in 
the case of the square lattice the S3Q phase has no net current flowing about any closed loops on the lattice but it has 
spin currents flowing around the elemental triangular plaquettes that exist in planes with Miller indices (1,1,1) (and 
equivalent symmetry planes) . The spin currents that flow correspond to a spin- flux of 7r/2 per triangular elemental 
plaquette (not 7r per plaquette as was the case in the square lattice). It is intriguing to note that Hasegawa et al. 
have pointed out that for a triangular lattice the optimal flux per plaquette in a U(l) flux phase is 7r/2 at (n) = 0.5 
[ 15 1 . Both the Siq and the S3Q states have been observed in 7-Mn alloys produced by doping with Fe, Ni, or Cu 
[24,H| and it would be of interest to see if spin currents can be detected in the S3Q phase of these materials. 

In conclusion, we have given numerical evidence that a spin-flux phase exists as a ground state of the Kondo lattice 
model with classical localized spins on a square lattice. This phase gives rise to a spin-flux of n for electrons circulating 
an elementary plaquette of the square lattice and manifests itself as a double-Q magnetic order in the classical spins. 
We have also proposed that a spin-flux phase may be stable on a FCC lattice. This phase manifests itself as a triple-Q 
magnetic order and gives rise to a spin- flux of ir/2 for electrons circulating the elementary triangular plaquettes that 
lie in the planes with Miller indices (1,1,1) (and equivalent symmetry planes). 
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